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We develope a theory of a fundamental effect of the interaction-induced decoherence of the 
electron wave function in a disordered metal. With the aid of the Keldysh technique and the path 
integral formalism we derive a formally exact equation of motion for the electron density matrix 
in the presence of interaction. We demonstrate that the effect of interaction of the electron with 
other electrons and lattice ions in a disordered metal is equivalent to that of an effective dissipative 
environment. Quantum noise of this environment causes quantum decoherence even at T = 0. 
Our analysis explicitely accounts for the Pauli principle which plays an important role for inelastic 
scattering processes but turns out not to affect quantum decoherence. Our results seriously challenge 
the existence of strong localization in low dimensional disordered metals. 



I. INTRODUCTION 



Recent experimentsEl attracted a great deal of attention to an old but fundamental question: how fast can a quantum 
particle loose information about its initial state in the presence of interaction? In other words, how fast can interaction 
destroy the quantum phase coherence? The answer on this question essentially depends on the type of interaction. 

■ It follows from general principles of quantum mechanics that quantum coherence of the wave function cannot be 
^vq , destroyed due to elastic interaction with an external potential. Another physical situation may take place if the 

quantum particle interacts with other (quantum) degrees of freedom which play the role of an effective environment. 
C — ■ ' In this case quantum dynamics of the particle cannot be described by the wave function but only by the density 
matrix. Various examples (to be discussed below) show that such interactions may lead to a complete destruction of 

■ quantum coherence. 

p3 ■ A general approach to the problem was formulated by Feynman and VernoroB who demonstrated that the effect 
of environment can be taken into account by means of averaging over its all possible quantum states. As a result 
the environment variables are integrated out and quantum dynamics of the particle can be described only in terms 
' of its own degrees of freedom. Within this approach interaction with the external environment is taken into account 
Q by means of the so-called influence functional which appears in the (effective) action for the particle as a result of 
O averaging over the bath variables. It is quite clear that specific properties of the environment are not important 
^ ^ unless they explicitely enter the expression for the influence functional. In other words, the particle does not "feel" 
•'H . the difference between physically different bathes provided thev,are described by the same influence functional. 

' These ideas were developed further by Caldeira and Leggetta who showed that the above arguments can be used 
^ . to describe quantum dynamics of dissipative systems and derived the rftective action for the case of linear Ohmic 
- - ' dissipation. The same type of analysis was also developed by SchmidB who formulated a quasiclassical Langevin 
equation approach describing real time dynamics of a quantum particle in the presence of dissipation and quantum 
noise. 

Although the above papers are dealing with the model of a bosonic environment it is obvious that the ideaad'tl 
can be applied to a fermionic bath as well. This was done e.g. by Ambegaokar, Eckern and Schorfl in the case of 
superconducting tunnel junctions and later by Schon and one of the present authorsQfl in a somewhat broader context 
of a metallic system with dissipation. Although the microscopic Hamiltonian describing electrons in a metafrcl is 
quite different from one used in the modda the final expressions for the influence functionals for various metallic 
systems and tunnel junctions obtained inETcl turn out to be similar or even completely equivalent to those considered 
in Refs. |4|j|. This equivalence is just an illustration of the property discussed above: the effect of physically different 
environments is indistinguishable provided they are described by the same influence functional. 

In Refs. ||-|| quantum dynamics of a certain collective variable of interest (the phase) was considered. This variable 
was extracted "from interaction" , after that electronic degrees of freedom were integrated out and the effective action 
for the phase was derived. One can also |Seneralizc this procedure and describe quantum dynamics of superconductors 
considering the phase as a quantum fieldo. In all these cases the collective variable is intimately linked to the electronic 
bath, quantum dynamics of the former does not exist without the latter at all. 
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In this paper we will analyze a somewhat different situation. Namely, we will study quantum dynamics of an electron 
propagating in a disordered metal and interacting with other electrons which play the role of effective environment. 
It is well known that quantpinJnterference of electrons scattered on impurities lead to quantum corrections to the 
classical Drude c0Hductivityll3O. These so-called weak localization corrections have been extensively discussed in the 
literature (see e.g.ll3ll3 for review). The magnitude of these corrections is known to be determined by the time within 
which electrons in a metal can be described by a phase coherent wave function. At times exceeding this so-called 
decoherence time r^p quantum coherence is destroyed, quantum interference is not anymore possible and therefore the 
classical diffusion picture is restored. 

The decoherence time in a disordered metal can be determined by various physical processea^^ch as electron- 
electron and electron-phonon interacticpa, electron scattering on magnetic impurities etc. (see e.g.cjiij). It was shown 
by Altshuler, Aronov and Khmelnitskiiiil that at not very low temperatures the effect of electron-electron interaction 
on tlie. .dccQherence time is equivalent to that of classical Nyquist noise in a disordered conductor. In this case one 

oc T^/f'^-^), (1) 

where d is the effective system dimension. This result demonstrates that decoherence effect of electron-electron 
interaction becomes weaker as the temperature is lowered. 

Down to which temperature does the above result remain correct? Or, more generally, does Tip increase with 
decreasing T at all temperatures thus going to infinity at T ^ 0? A positive answer on the latter question would 
mean that at T = the electron in a disordered metal can be described by the phase coherent wave function even in 
the presence of Coulomb interaction with other electrons. This scenario faces several important problems: 

(i) Even if one assumes that electrons in a metal are isolated from the external environment it is hard to believe 
that the iV-electron wave function ^I^jv {N = 10" is the total number of electrons in a conductor, a ^ 1) can be 
factorized into the product of single electron wave functions. This would imply that Coulomb interaction is effectively 
renormalized to zero at T = 0. We see no reason to expect that. 

(ii) In any real physical situation electrons in a metal are coupled to the external environment, e.g. via external 
leads. Electrons can enter and leave a piece of metal, interact with electrons in the leads etc. Due to this reason the 
whole system of N electrons in a metal (and not only one electron!) should be treated as open and thus cannot be 
described by the wave function ^jv even at T = 0. 

(iii) In the case of Caldeira-Leggett type of models it is well known that after a sufficienly long time even at T p. 
quantum coherence is destroyed due to interaction with the bath no matter how weak this interaction is (see e.g.llll). 
For strong enpugh interaction the particle exhibits practically no quantum features and can even get localized in a 
potential wcllEa-E3. Since the Caldeira-Leggett models are qualitatively similar to the problem in question (in both 
cases one quantum degree of freedom interacts with many others playing the role of a bath) it is quite natural to 
expect also the decoherence time for interacting electrons in a disordered metal to remain finite even at T = 0. 

(iv) Saturation of the temperature dependeujce jEor at low T was clearly observed in experiments^ as well as in 
numerous earlier experimental studies (see e.gBHcJ). Even a qualitative interpretation of these results in terms of 
eq. (0) is not possible. There exist sufficiently strong experimental reasons to believe that this saturation is not due 
to heating of the sample or the effect of magnetic impurities. 

All the above arguments strongly suggest that the decoherence time should remain finite down to T = 0. 
Although the experimental evidence for the low temperature saturation of is available for quitp-|a long time the 
physical reason for this effect was not convincingly explained in the literature. An earlier attempted to promde this 
explanation in terms of zero-point fluctuations of impurities has failedEj. Only very recently it was sugg estedtlO that 
quantum decoherence in a disordered conductor is caused by the mechanism qualitatively similar to one existing in 
the Caldeira-Leggett model, namely the interaction of a particle with an external quantum bath. In the present case 
an electron propagating in a disordered metal interacts with a fluctuating electric fleld produced by other electrons 
which play the role of an effective dissipative environment. Due to randomness of such fluctuations the quantum 
coherence is lost within a finite time t^. 

In fact, this mechanism is not new. Exactly the same effect was considered in Ref. ^ where only the low frequency 
(classical) fluctuations of the environment uj < T were taken into account. This approximation is sufficient at high 
temperatures and yields the result (|l]). However, as the temperature is lowered there remain less and less classical 
modes in the fluctuating environment. Therefore it is reasonable to expect that below a certain temperature Tq 
quantum fluctuations of the environment W|iili.j|U_ji_.T will take over aud determine the decoherence time r^. The 
generalization of the corresponding analysisEjEjO'Ej is straightforwardEj and allows to |evaluate both the classical- 
to-quantum crossoveij-temperature Tq and the decoherence time at low T. The resultfcH agree quite well with the 
experimental findinggj. 

Several fundamentally important questions should still be answered, however. Perhaps the main one is about the 
role of the Pauli principle in our description of the electron-electron interaction. At sufficiently high temperatures 
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there are always enough empty states into which electrons can get scattered provided the energy transfer due to 
interaction does not exceed T. In this case the Pauli principle does not play a very important role. On the other 
hand, at T ^ in equilibrium all quantum states below the Fermi energy are occupied. Thus scattering into any of 
these states is forbidden due to the Pauli principle and the electron energy cannot decrease. At T = it also cannot 
increase due to interaction with other electrons because they are also in the ground state. Thus at T = inelastic 
electron-electron scattering is forbidden due to the Pauli principle. If so, how can one expect quantum decoherence 
to occur as a result of the electron-electron interaction at T = 0? One might think that due to the Pauli principle the 
electron is simply insensitive to zero-point fluctuations of other electrons, and quantum decoherence does not take 
place. 

Another though somewhat related question is about the effectpof phonons. It is quite obvious that at a finite T 
decoherence can also be caused by electron-phonon interactioi£3ll3. However at T ^ there are no real phonons 
in the system and therefore no quantum decoherence due to phonons should be expected. On the other hand, if we 
admit that zero-point fluctuations of electrons lead to decoherence, zero point fluctuations of the lattice ions should 
have qualitatively the same effect. Is it possible to match this conclusion e.g. with the well known fact that electrons 
do not scatter on zero-point lattice fluctuations, and therefore the phonon contribution to the conductance of a metal 
vanishes at T = 0? 

On a more formal level it is sometimes conjectured that the above arguments correspond to the exact cancellation 
of Keldysh diagrams for Tip in the limit T — > 0. If this is indeed the case the inverse decoherence time 1/T(^ would be 
exactly equal to zero at T = and (at least qualitatively) would match with the result (Q) . 

The main goal of this paper is to demonstrate that this is not the case. 

The key point is that all the above arguments apply as far as the electron inelastic relaxation time is concerned. 
The latter can be found by means of a quantum kinetic approach which will be fully reproduced| W |i tl j i | in the framework of 
our analysis. However, as it has been pointed out in Ref. ij and discussed later in other worksc§ll30, the decoherence 
time is entirely different from and, in contrast to the latter, cannot be determined with the aid of any kind of 
kinetic analysis. Roughly speaking, the decoherence time controls the decay of off-diagonal elements of the electron 
density matrix and (even at high temperatures) has little to do with inelastic scattering of electrons. Therefore 
the quasiclassical kinetic treatment is insufficient and a fully quantum mechanical analysis should be elaborated to 
determine the time t^. 

The paper is organized as follows. In Section 2 we make use of the general formalism of the Green-Keldysh functions 
and derive a formally exact equation of motion for the electron density matrix in the presence of Coulomb interaction. 
This equation explicitely accounts for the Pauli principle and allows for a clear understanding of its role in the 
process of electro-electron interaction in a metal. In Section 3 the effective action (or the influence functional) for the 
fluctuating scalar potential in a metallic conductor is derived. We also demonstrate that in equilibrium this influence 
functional satisfies the fluctuation-dissip|atjpi|i theoremEJ and establish the relation with the real time effective action 
derived in the Caldeira-Leggett modelgQ'ErElnJ. In Section 4 with the aid of these general results we will derive the real 
time effective action for the electron propagating in a metal and determine the decoherence time r^p and the weak 
localization correction to conductivity at low temperatures. Our formalism naturally includes both electron-electron 
and electron-phonon interactions and allows to establish the corresponding contributions to from each of these 
processes. In Section 5 we derive the quasiclassical kinetic equation and demonstrate the relation of our analysis to the 
standard kinetic approach which allows to evaluate the inelastic scattering time r^. We also derive the quasiclassical 
Langevin equation which under certain conditions can be used to describe propagation of electrons in a disordered 
metal. Discussion of the results is given in Section 6. We discuss the relation of our results with previous works as 
well as possible consequences of the low temperature saturation of for the existing picture of strong localization 
in low dimensional conductors. We also briefly compare our theoretical predictions with available experimental data. 
Some details of our calculation are presented in Appendices. 



II. DENSITY MATRIX 



We will consider a standard Hamiltonian describing electrons in a disordered metal 



Hei = Hq + Hint 



(2) 



where 




(3) 
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Here is the chemical potential, U{r) accounts for a random potential due to nonmagnetic impurities, and v(r) — l/\r\ 
represents the Coulomb interaction between electrons. 

Let us define the generating functional for the electron Green-Keldysh functions in terms of the path integral over 
the Grassman fields V' and ■0 

^ JVVjVi^JVi^ exp {iSeff [yi, 0, F] + z dt J drm, r)7?(t, r) + ^{t, r)rf {t, r]) 

where Seff is the effective action 

Seff[i>, V, = y^/* (/ ^^[*^(^' ^Wi*^ r) - e^{t, r)^{t, r)V{t, r)] ~ H^. ^[ 



+ 



dtjdrj dr'V{t,r)v-^{r - r')V{t,r'), (6) 



where v~^{r — r') = — V^/47r. Integration over time_f in (o) goes along the Keldysh contour K which runs in 
the forward and then in the backward time directionsEJ. In (|^,^) we performed a standard Hubbard-Stratonovich 
transformation introducing the path integral over a scalar potential field V in order to decouple the i/i^-interaction 
in (^) . The electron Green-Keldysh function G can be determined from (|^) by taking the derivatives with respect to 
the source fields i] and rj* : 



S 

Si] 

Making use of (||, 0) and the definition of the Green-Keldysh function for an electron interacting with the field V 



G{t,r;t',r') = z^-^^-^-^ J[r;, ,y*]|,=^.=o. (7) 



it is easy to prove the identity 



- ^ / Gy e'SlVuV2] 

j VV1VV2 eiS[Viy2] ' 

where 

iS[Vi,V2] = 2Ti\nGy^ +1 J dt' J dr ^^^^^ g^^^^"^ ' 



The factor 2 in front of the trace comes from the summation over a spin index. In 10) we explicitely defined 
the fields Vi{t) and V2{t) equal to V{t) respectively on the forward and backward parts on the Keldysh contour K. 
Analogously G and Gy = G[Vi, V2] are the 2x2 matrices in the Keldysh space: 

(11) 
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The matrix function Gy obeys the equation 

._d_ 



z J- - iJo(ri) + eV{ti,ri) ] Gy = 5{ti - t2)5{ri - rj); (12) 



where 
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Note that the function Gy is to some extent similar to the Green-Keldysh function of an electron in an external field. 
However there exists an important difference: in our case the electron interacts with a fluctuating (quantum) field V . 
Formally this implies that the fields V{t, r) on two parts of the Keldysh contour differ Vi{t, r) ^ V2{t, r) while for the 
external field one always has Vi(t, r) = VAt^r). 

The general solution of the equation (n2h can be expressed in the form 



Gv{ti,t2) = -iUv{ti,t2)[e{h ^ t2)a - e{t2 - ti)b + fv{t2)], 



Here we defined 











I 


j ' 







Uv(ti, t2) is the matrix evolution operator 

Uv{ti,t2) = 

which consists of the scalar evolution operators 
ui.2{ti,t2) = Texp 



Ul{ti,t2) 

U2{ti,t2) 



(14) 



(15) 



(16) 



i J dt'(Ho-eVi^2{t')) 
ti 

r *^ 

J Vr{t')exp i J dt' 



r{t\)—ri 



mr 



U{r)+eVi,2{t\r) 



(17) 



T is the time ordering operator. In eq. (|lj) and below we always imply integration over the internal coordinate 
variables in the product of operators, whereas integration over time is written explicitely. For the sake of brevity we 
also do not indicate the coordinate dependence in ( p^ ) and many subsequent expressions. This dependence can be 
trivially restored if needed. 

Note that eq. (|l^) is completely equivalent to the standard representation of the Green-Keldysh matrix which 
elements can be expressed in terms of retarded, advanced and Keldysh Green functions. The representation (|lj) 
defines a general solution of the linear differential equation (|l^): the term Uv{ti,t2)fv{t2) with an arbitrary matrix 
operator fv{t2) represents a general solution of the homogeneous equation, while the terms with 0- functions give 
a particular solution of the inhomogencous equation. The operator function fv{t2) m (Q) is fixed by the Dyson 
equation 



Gv{tut2) = Go{tiM)- / dt'GQ{h,t')eV{t')Gv{t'M)- 



(18) 



The matrix Go is the electron Green-Keldysh function without the field. This function is defined by eqs. ( p^| - |l7| ) with 
Vi,2{t, r) = and /o(t2) has the form 



/o(i2) 



1 - Po(i2) Poit2) 



(19) 



where po{t) = e *^f'p(0)e is the electron density matrix for Vi,2 = at a time t. 

The equation (|l|) can be solved perturbatively in eV. Combining this solution with eq. (^ one reproduces the 
standard Keldysh diagrams. This way of treating the problem is quite complicated in general and becomes particularly 
nontransparent in the interesting limit of low temperatures. 

We will proceed differently. 

It is well known that the 1,2-component of the Green-Keldysh matrix G is directly related to the exact electron 
density matrix 



p{t;r,r') = -iGi2{t,t;r,r'), 



(20) 
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which contains all necessary information about the system dynamics in the presence of interaction. Analogously 
one can define the "density matrix" pv{t) = fi,2{t) related to the 1,2-component of the matrix Gv by the equation 
equivalent to (^0|). Our strategy is as follows. First we will derive the exact equation of motion for the density matrix 
pv{i) which turns out to have a very simple and transparent form. Already at this stage we will clarify the role of the 
fluctuating fields Vi^2 and the Pauli principle in our problem. Then we will evaluate the influence functional ^[V^i, V2] 
( |lo| ) and find the density matrix p from the equation 

p(t;r,r') - (/5y(t;r,r')y,,y,, (21) 

where the average over fields Vi and V2 is defined in (^). 

The derivation of the equation for the density matrix pv{i) is straightforward. Let us perform the time integration 
in the last term of eq. (|l8|) . Integrating by parts and making use of eq. (jlj) after a simple algebra (see Appendix A) 
we obtain 

Go{ti , t)G{t, <2) - Go(ii , 0)G(0, ta) - 0. (22) 

Substituting the representation ( p^ into ( p2| ) we arrive at the matrix equation which relates the matrix fv{i), the 
evolution operator Uy{t) and the initial density matrix Pq{Q) defined for Vi,2 = 0. With the aid of this equation one 
determines the 1,2-componcnt of the matrix fv{t) and thus the density matrix pv(t). The details of this calculation 
are presented in Appendix A. As a result we find 

[1 - PQit){u2{t, 0)ui(0, t) - l)]pv{t) = po{t). (23) 

One can also rewrite this result in the form of the differential equation describing the time evolution of the density 
matrix: 

t^^[Ho,pv]-{l-pv)eVipv+PveV2il-pv), py(0) = po(0). (24) 

The equation ( |2^ ) is the main result of this Section. We would like to emphasize that our derivation was performed 
without any approximation, i.e. the result ( p^ ) is exact. It contains all information about the system dynamics hidden 
in the four components of the Green-Kcldysh matrix. In the absence of the fluctuating field Vi_2 = eq. ( p^ ) reduces 
to the standard equation for the electron density matrix with the Hamiltonian Hq. In the presence of the field ^1^2 
the equation (|2j) exactly accounts for the Pauli principle. This is obvious from our derivation which automatically 
takes care about the Fermi statistics through the integration over the Grassman fields tp and "tp. This is also quite 
clear from the form of the last two terms in the right hand side of eq. (p^. In Section 5 we will demonstrate that 
within the quantum kinetic analysis these terms are responsible for the standard in- and out-scattering terms in the 
collision integral. 

In order to understand the role of the Pauli principle let us rewrite the equation (|2j) in the form 

= [^0 - eV+,py] - (1 - pv)^pv - Pv^il - Pv), (25) 

where we defined = {Vi +V2)/2 and — Vi — ¥2- It is quite obvious from ( p5| ) that the field V^{t, r) plays the 
same role as an external field. All electrons move collectively in this field, its presence is equivalent to local fluctuations 
of the Fermi energy /x — *■ /i -I- eV~^{t, r). The Pauli principle does not play any role here. Below we will demonstrate 
that quantum fluctuations of the field y+ are responsible for the low temperature saturation of the decoherence time 
and the weak localization correction to conductivity in disordered metals. 

The field V~ is, on the contrary, very sensitive to the Pauli principle. It will be shown below that this field is 
responsible for damping due to radiation of an electron which moves in a metal. Corresponding energy losses can be 
only due to electron transitions into lower energy states. At T = in equilibrium all such states are already occupied 
by other electrons, therefore such processes are forbidden and the electron energy remains unchanged due to the Pauli 
principle. We will demonstrate, however, that these processes are irrelevant for the decay of the off-diagonal elements 
of the electron density matrix and therefore the Pauli principle can hardly affect quantum decoherence in a disordered 
metal even at T = 0. 

III. INFLUENCE FUNCTIONAL FOR THE FIELD 

Let us now derive the expression for the influence functional (effective action) S[Vi, V2] for the fleld V. A formally 
exact action S[Vi, V2] obtained by integration over all electron degrees of freedom is given by eq. (|lO|). Let us expand 
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this expression up to the second order in V1.2. The first order order terms of this expansion vanish because the Green 
function Gq corresponds to a zero current and zero charge density state of the system. In the second order we obtain 

iSS^^^ = -TiiGoeVGoeV) 



= -e^tr 



G11V+G11V+ - G12V+G21V+ ~ G21V+G12V+ + G22V+G22V+ 



GiiV+GiiV~ - Gi2V+G2iV~ + G21V+G12V- - G22V+G22V^ 
^GnV-GnV- + ^Gi2V-G2iV- + ^G2iV-Gi2V- + iCaay^GaaV^" 



(26) 



Here the subscript for the Green functions is omitted for the sake of simphcity, all the Green functions here and 
below in this Section are defined for V1.2 = 0. The expression ( p6| ) can be simplified with the aid of the identities 

Gil = G12 + G^ = G21 + G^, G22 — G12 — G^ — G21 — G^ , (27) 
which allow to exclude the functions Gn and G22 from the action (|2^). The terms containing V^V^ are reduced to 

tr{G^V+G"V+ + G'^V+G'^V+). (28) 

Making use of the expressions 

G"^ -ie{h-t2)uoih,t2), G^ = ie{t2-h)uoih,t2), (29) 



and writing the traces ( P8[ ) in the time-space representation, we immediately observe that the product 9{ti —^2)^(^2 — 
ti) = appears under the integral, and the whole combination ( p8| ) vanish. In other words, the terms V~^V~^ give no 
contribution to the action (p6|). The remaining terms in combination with the last two terms in (|l^) describing the 
free field action yield 



tS[Vi,V2] 



I J dt' I dr 





47r 



t t 
dti J dt2 I dridr2 x 





iV-{ti,ri)x{ti,t2,ri,r2)V+{t2,r2) + V {ti,ri)T]{ti,t2,ri,r2)V (i2,r2) 



(30) 



where 



X{ti,t2,ri,r2) = -ie'^[G'^iti,t2,ri,r2){G2i{t2,ti,r2,ri) + Gi2(i2, ^i, ^2, ri)) ^ 

{Gi2{ti,t2,ri,r2) + G2i{ti,t2,ri,r2))G^{t2,ti,r2,ri) 



(31) 



V{ti,t2,ri,r2) = — [Gi2{ti,t2,ri,r2)G2i{t2,ti,r2,ri) + G2i{ti,t2,ri,r2)Gi2{t2,ti,r2,ri)] 



(32) 



The expressions (|30|-^2|) define the influence functional for the field V in terms of the Green-Keldysh functions 
for the (in general nonequilibrium) electron subsystem. It is easy to check that the above expressions satisfy the 
requirement of causality: in the V~ {ti)V'^ (t^) terms the time ti is always larger than t2 due to analytic properties 
of retarded and advanced Green functions (29). It is also straightforward to demonstrate (see Appendix B) that in 
thermodynamic equilibrium the kernels ri{uj) and Imx(a;) (Fourier transformed with respect to the time difference 
ti — ^2) satisfy the equation 



r](uj,ri,r2) 



^ coth ( ^ j Imx ( w , r 1 , r 2 ) . 



(33) 



The latter equation is just the fluctuation-dissipation thcorem0. 

Finally, for a homogeneous system one can also perform the Fourier transformation with respect to ri — r2 and find 



iS[Vi,V2]^i I -j^:;^^ {~^,-k) — ^ — v+{uj,k) 



dujd^k^^_, , ^ fc2lme(tj, fc) , / w \ , ,x 
^ 7^^coth(— jF-(^,fc); 



(34) 
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where e{uj, k) is the dielectric susceptibiHty of the system. For a homogeneous electron gas it is given by the standard 
RPA formula: 



e{u, fc) = 1 + 4- / d^p (35) 

7r2fc2 J ^ u - S^p+k + S.P + 

Eq. (|3^) includes only the electron contribution to the susceptibility. In general the effect of ions should also be 
accounted for. Here we will describe this effect within a very simple approximation which is, however, sufficient for 
our analysis. Namely, the ion contribution to the susceptibility will be taken in the form 5ti = — Wpj/a;^, where ojpi 
is the ion plasma frequency. Then the phonon spectrum is determined by the equation e(a;, k) = 0. In the long wave 
limit this approximation works sufficiently well for longitudinal phonons which mainly interact with electrons. The 
effect of transverse phonons cannot be described within this simple model. But such phonons are weakly coupled to 
the elecrons anyway, and therefore their effect can be safely ignored. If needed, further generalizations of this simple 
model can be also incorporated into our analysis. i-. 

In the relevant case of a disordered metal a direct calculation of polarization bubbles (see e.g.El) yields 

<^,fc) = l+ .^7n., -4- (36) 
—lu + Dk^ uj'' 

Here tr — 2e^Ni^D is the classical Drude conductivity, Nq is the metallic density of states and D ~ vpl/3 is the 
diffusion coefficient. 



The expression (36) is valid for wave vectors smaller than the inverse elastic mean free path k '^l/l and for small 



frequencies w < I/te, where Te = l/vp is the elastic mean free time. Note, that if one neglects the effect of phonons 
and considers only nearly uniform in space (fc « 0) fluctuations of the field V |0|Be| .immediately observes that eqs. 
(^,^6|) exactly coincide with the real time version of the Caldeira-Leggett actionQflu^El in this limit. For a disordered 
metal this action was derived by means of the quasiclassical Eilenberger equations in Refs. ^|J|. Our analysis (see 
alsoB) reproduces these results and generalizes them to the case of nonuniform fluctuations of the field V . 

Taking into account only uniform fiuctuations of the electric field one can also derive the Caldeira-Leggett action 
expressed in terms of the electron coordinate only. In the quasi-one-dimcnsional situation one should simply write 
down the electron action on the Keldysh contour, take into account the potential ener gy —Ex of the electron with 
the coordinate x{t) in the fiuctuating electric field E{t) and add the action for the field (34 3^) (with the last term in 



( pq ) being dropped). After identification E{uj) = —ikV{LJ,k) — > —V{lj)/L (L is the sample length) and integration 
over the fluctuating field V one arrives at the Caldeira-Leggett action for the electron coordinate x{t). In this case 
the effective viscosity in the Caldeira-Leggett influence functional is proportional to 1 /a (in contrast to the effective 
viscosity for the field V which is proportional to a). The whole procedure is completely analogous to that discussed 
in details in Ref. ^ where we considered the real time effective action for a dissipative system characterized by two 
collective degrees of freedom (the phase and the charge). Integrating over the charge variable one arrives at the 
Caldeira-Leggett action in the "phase only" representation. The same can be done here if we use a formal analogy of 
x{t) with the phase and V/iuj with the charge (as defined in Ref. p9| ). 

For our present purposes it is not sufficient to restrict ourselves to uniform fluctuations of the collective coordinate 
V of the electron environment. We will see that fluctuations with nonzero k play an important role and should be 
taken into account in the quantitative analysis. The corresponding effective action will be derived in the next section. 
However the main message is clear already from the simple example considered above: in a disordered metal the effect 
of Coulomb interaction of the electron with other electrons is equivalent to that of an effective dissipative environment 
with the correspondent effective viscosity governed by the Drude conductivity a. 

IV. DECOHERENCE TIME AND CONDUCTANCE 

A. Conductance and electron effective action 

In order to evaluate the system conductance we will determine the single particle density matrix p in the presence 
of an external electric potential Vx{r) applied to the metal. Generalization of the results obtained in the previous 
sections to the case 14 7^ is straightforward. The density matrix is determined by the equation 

Pit) - {pv{t,Vx)))v,,v, jp^^p^^g.s.jv.^y.] ' (37) 

where the effective action 5'y^[Vi, V2] reads: 
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I J dt' I dr 





8tt 



(38) 



where the subscript V + Vx indicates the shift of the fields V1.2 ^ Vi,2 + Vx- Here we use the formahy exact expression 
for the effective action, the approximation (Q) will be introduced after the expansion in will be carried out. 

The density matrix pvit, Vx) obeys the equation ( ^ ) with + Vx- Assuming the field Vx sufficiently small 

one can linearize the equation for pvit, Vx) = pv{t) + Spv{t, Vx) and get 



.d5p\, 



HiSpv - SpvH2 - [eVx,pv]; 



(39) 



where 



Hi^Ho- eV^ 



(1 - 2pv)eV^ 



H2 = Ho- eV+ + -eV-{l - 2p 



(40) 



The formal solution of the equation (39) can be easily found: 

t 

dpv{t)^i j dt'Ui{t,t')[eVx.pv{t')]U2{t\t), 


where 



(41) 



C/i,2(ii,i2) = Texp 



dt'Hi^2{t') 



(42) 



The operators Hi 2 (40) are nonlocal (since they contain the density matrix), therefore the path integral representation 
for the evolution operators ( ^ ) contains an additional integration over momentum. The operators (^) can be written 
in the form: 



i?i(p,r) - 1^ + U{r) - eV+{t,r) 2n(i/o(p, r))] ey-(i, r), 



2m 

^2 



H2{p, r) = ^ + U{r) - eV+it, r) + \eV-{t, r) [l ~ 2n(i/o(p, r))] ., 



(43) 



where n.(^) = l/[exp(^/r) + 1] is the Fermi function. In deriving (^) from ( ^o|) we set pv to be an equilibrium 
density matrix. We also neglected the effect of Coulomb interaction in the expression for pv{t)- This approximation 
is justified as long as Coulomb interaction is sufficiently weak. Note that the same approximation for py(t) can be 
used in eq. (pl|). 

The evolution operators (^2|) acquire the form: 



r(t2)=rf r t2 

Ui,2{tiM;rf,ri)= j Vr{t') J Vp{t') exp I J dt'{pr- Hi^2{p-r)) 



r{ti)—ri 



(44) 



For the sake of generality we note that in the presence of interaction there exists an additional (linear in the field) 
correction to the density matrix. In order to see that let us expand the action (^8|) to the first order in Vx'. 



iSS[Vx] = -2etr((Gii - Gza)^.). 



(45) 



This correction to the action gives an additional contribution to the density matrix (|37|). Expressing the functions 
Gil and G22 in terms of the density matrix pv and the evolution operators ui^2 and combining this correction to the 
density matrix with one defined in (|l]) we find 



Sp(t) = {Spv)v+,v-- + {Spint)v+,v- 



(46) 
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where 5pv is given by (^Tj) and ^pint has the form 

t 

5pi„t = -27;py(t) [ dt'iY{u^{t,t')[eV^{t'),pv{t')]u2{t\t)). (47) 



In the hmit of weak interaction between electrons the averaging in (^) may be performed with the approximate 
action S[Vi,V2] (|3|). 

It is easy to observe that the second term in (^6|) is small in the limit of weak interaction and vanishes completely 
if interaction is neglected. The weak localization correction is described by the first term in eq. ( p6| ) which will be 
only considered further below. 

Making use of a standard definition of the current density j: 

j{t,r)^- {Vr,5p{t,ri,r2)-Vr,5p{t,r^,r2))\,^^,^^,, (48) 

combining it with the above equations and assuming the external electric field to be constant in space and time, 
Vx ~ —Er, we arrive at the expression for the system conductance 

^ t 

CT=^ J dt' J dradri2 {"^nf - '^r2f)\^^^^^^^ J{t,t'-rif,r2f;rii,r2i){rii - r2i)po{ru,r2i). (49) 

— OO 

Here we have shifted the initial time to — cxd. The function J is the kernel of the operator 

j^J2u\v){v\u+ 

V 

where the sum runs over all possible states of the electromagnetic environment. This function can be expressed in 
terms of the path integral 

J{t,t';rif,r2f]rii,r2i)^ J Vri j 'Dr2 J 'DpiVp2 x 

ri(t')=rii r2{t')=r2i 

JSo[ri,pi]-iSo[r2,P2]+i J^, dt" J dr{f-V++f+V-) 



(50) 



v+,v- 



Here the action So[x,p] has the form 



So[r,p] = J dt" (^pr - ^ - Uir)^ , (51) 
t' 

and the "charge densities" /~, are defined by the equations: 
nt, r) = eSir - r,{t)) - eS{r - r2{t)), 

f+{t,r) ^l(e[l- 2nip,{t),riit))]S{r - r,{t)) + e[l - 2nip2{t),r2it))]S{r - r2{t))^ . (52) 

Averaging over , V~ in eq. ( ^ ) amounts to calculating Gaussian path integrals with the action (^J) and can be 
easily performed. We obtain 

ri{t)=rif r2{t)=r2f 



J{t,t']rif,r2f;ru,r2i)= j Vri J 'Dr2 J 'Dpi'Dp2 x 



ri{t')=rii r2{t')=r2i 

X exp{«5o[ri,pi] - iSn[r2,P2] - iSR[ri,pi,r2,P2] - S'/[ri,r2]}; (53) 

where 
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t 

„2 



5fl[ri,pi,r2,P2] - Y J dh J dt2{Riti-t2,ri{ti)-ri{t2))[l-2n{pi{t2),ri{t2))]- 
t' t' 

-R{h ~ t2,r2{ti) - r2(t2)) [1 - 2n{p2{t2), r2(t2))] 
+R{ti - i2, r-i(ii) - r2(t2)) [1 " 2n{p2{t2), r2{t2))] ~ 

-R{ti - t2,r2{ti) - ri(t2))[l - 2n(pi(t2),ri(t2))]}; (54) 



and 



t t 



Si[ri,r2] = Y J dh J dt2!^I{h ^ t2,ri{h) ~ ri{t2)) + Hh - t2,r2{ti) - r2{t2)) - 
t' t' 

- t2,ri{h) - r2(i2)) - I{h ~ t2,r2{ti) - ri(t2))|. (55) 
At the scales > Z the functions R and / are defined by the equations 



UJ 



dujd^k / -An 

If necessary, more general expressions for R and / for the whole range \r\ > l/pF can be easily derived. We will avoid 



Iit,r) = I :^lnr J coth )e— . (57) 



doing this here because the approximation (57) is already sufficient for our present purposes. 

Note that the expression in the exponent of eq. (^ ) defines the real time effective action of the electron propagating 
in a disordered metal and interacting with other electrons. The first two terms represent the electron action Sq ( |51| ) 
on two branches of the Keldysh contour while the last two terms Sn and Si determine the influence functional which 
comes from the effective electron (and/or phonon) environment. As can be seen from eqs. (p^^7|) this influence 
functional is in general not identical to one derived in the Caldeira-Leggett model. However on a qualitative level 
the similarity is obvious: in both models the influence functionals describe the effect of a certain effective dissipative 
environment. 

Let us now analyze the expression ( p9| ) and establish connection with other conductivity calculations. It is convenient 
to introduce the Wigner function n{p,r) — J (ir~e~'f"' p{r + r~/2,r — r~/2) instead of the density matrix. For 
homogeneous systems it does not depend on r at the scales exceeding the mean free path. With this in mind we 
obtain 



where W{t^ t') is the evolution operator for the Wigner function: 

n{t, X, p) ^ W(t, t')n(t\ X, p). (59) 

The kernel of this operator and that of J (^3|) are related by means of the Fourier transformation with respect to 
'''14,1/ — '''24, 2/- In equilibrium one has ri(p) = l/(exp(^/T) + 1), and therefore at small T one finds 

and arrives at the standard result (cf. e.g.El) 

^ t 

^^2e_A^ j dt'(v(t)v(t'))- (60) 
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B. Weak localization correction and decoherence time 



Let us analyse the structure of the function J (|5^) in the same spirit as it has been done in Ref. y_3|. In the zero 
order approximation one can neglect the terms Sn and Si describing the effect of Coulomb interaction. Then in the 
quasiclassical limit ppl 3> 1 the path integral (^3|) is dominated by the saddle point trajectories for the action Sa 
which are just classical paths determined by the Hamilton equations 



P 



dHo{p,r) 
dr 



r = 



dHoip,r) 
dp 



(61) 



with obvious boundary conditions ri{t') = ru, r{t) = rif for the action S'o[t'i,Pi] and r2{t') — r2i, i'2{t) — r2f 
for the action 5*0 [^2,^2]- Substituting these saddle point trajectories into (|53| ) and integrating out small fluctuations 
around them one finds 



J{t, t'\ rif, r2f;rii, r2i) X] "^''i X! t'] ^i/, Vu) ~ iSo{t, t'; r2f, r2i)) 



(62) 



where the actions 'S'o(i, t'; Ti_2/, ^"1,21) are taken on the classical pathes Ti.2(i) and 




det 



\ drifdru 



(63) 



The value A*^ is defined analogously. 

Since in a random potential U{r) there is in general no correlation between different classical paths ri{t) and r2(t) 
these pathes give no contribution to the double sum (62): the difference of two actions in the exponent of (62) may have 
an arbitrary value and the result averages out after summation. Thus only the paths for which iS'o[ti,pi]~~ Sq Xc,p2] 
provide a nonvanishing contribution to (|6^). Two different classes of such paths can be distinguished (see e.g.E3): 

i) The two classical paths are almost the same: ri{t") ~ r2{t"), Pi{t") ~ P2{t") (see Fig. la). For such pairs we 
obviously have ru ~ r2i and rif r2/. In other words, in the path integral (53) one integrates only over trajectories 
with \ri{t") — r2(i")| ^ 1/pf- Physically this corresponds to the picture of electrons propagating as nearly classical 
particles which can be described by the diagonal elements of the density matrix only. In the diffusive limit these paths 
give rize to diffusons (see e.g.llHl) and yield the standard Drude conductance. 

it) The pairs of time reversed paths. In this case ru ~ r2/, r\j ~ r2i (Fig. lb). In the path integral (53) the 
trajectories r\ and r2 are related as r2{t") ~ ri(t + t' — t") and P2{t") — ~Pi{t + t' —t"). In other words, in (53) one 
integrates over paths with \r2{t") — ri{t + t' — t")\ < 1/pf, however the difference \ri{t") — r2{t")\ maybe arbitrarily 
large in this case. In a disordered metal these paths essentially determine the dynamics of off-diagonal elements 
of the electron density matrix. They correspond to Cooperons suad-give rize to the weak localization correction to 
conductivity. This correction is defined by the following equationllnij 



5(7 = - 



'^''^ J dtW{t) = -'^ J dt'W{t~t') 



2e^D 



(64) 



Here we changed the parameter of integration in order to make the relation with the equation (^9|) more transparent. 
The quantity W{t) represents the effective probability for the diffusive path to return to the same point after the 
time t. Note that W{t) contains the contribution from time reversed paths and therefore differs from a classical 
probability. However, in the absence of any kind of interaction which breaks the time reversal symmetry this value 
coincides with the classical return probability and is given by the formula Wo{t) = {4:nDt)-'^/^a-'^^''^\ where d is the 
system dimension and a is the transverse sample size (the film thickness for d = 2 and the square root of the wire 
cross section for d = 1). 

The weak localization correction (^) diverges for d < 2. This divergence can be cured by introducing the upper 
limit cutoff for the integral (|6^) at a certain time . This time is usually reffered to as decoherence time. As we have 
already discussed, in a disordered metal the time is determined by electron-i eLect ron, electron-phonon and other 
types of interaction which may destroy quantum coherence. From (p3) one findal3't3: 



V3e 
27r3/2;' 

■^In 



d = 3, 
d=2 
d=l. 



(65) 
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a b 




c 

FIG. 1. Different quasiclassical trajectories: a) ri and r2 are close to each otfier; b) time reversed paths; c) time reversed 
paths returning to the vicinity of the initial point, these paths determine the weak localization correction to conductivity 



Here and below — aa?'^'^ is the Drude conductance of a d-dimensional sample. 

To evaluate T^p we first note that the functions R and / (^) change slowly at distances of the order of the Fermi 
wavelength 1/pf. Therefore we may put r\(t") — r{t"), r2{t") — rt{t") = r{t + t' — t"). Here r(t") is a classical 
trajectory with the initial point r{t') = (its position can be chosen arbitrarily) and the final point \r{t)\ < I (Fig. Ic). 
In other words, we consider trajectories which return to the vicinity of the initial point. Then the return probability 
W{t) can be written in the form 

W{t ~ t') ~ Wo{t - t') ^e-*^«[*'*'''-'P''-*'P'l-^^[*'*'''-''-*lJ> , (66) 

where the average is taken over all diffusive paths returning to the initial point. The average in (|6^) decays exponen- 
tially in time, therefore we may define as follows: 

-- - ' -iSR[t,t':r,p:rt,Pt]-Si[t,t'-,r,rt] 



e ^'p = { e 



_ ^{-^Snlt.t'■.r,p■.r„p,]-S,lt,t'■r.r,])^^ _ _^ ^g^^ 

where the average is again taken over all classical paths returning to the same point at the time t. 

Let us consider the term Sr (|5^). Making use of the obvious relations R{t, —r) = R{t, r), n{—p, r) — n{p, r) after 
a trivial algebra we get: 

2 ' * 

SR[t,t';r,p;rt,pt] = ^ J dh J dt2l^2R{ti - t2,r{ti) - r{t2))[n{p{ti),r{h)) - n{p{t2),r{t2))] + 
t' t' 

+R{h +t2~t~ t',r{ti) - r{t2)) [1 - 2n{p(t2),r{t2))\ - 

-Riti +t2-t~t', r{t + t'^ h) - r{t + t' - t2)) [1 - 2n{p{t + t' - ia), r{t + t' ~ tz))] | (68) 

It is clear that the difference of the occupation numbers n in the first term of eq. ( |68|) vanish: n depends only 
on the energy and not on time because the energy is conserved along the classical path (|6l|). Thus the first term in 
the integral (^) is identically zero already before averaging over diffusive paths. The difference of two other terms 
vanishes after averaging over such paths. It can be easily seen if we again bear in mind that the occupation numbers 
do not depend on time and make use of the fact that the average {R{t3, r{ti) — r{t2)))r is a function of the time 
and the absolute value of the time difference |ti — t2\. This implies that after averaging the two last terms in eq. ( pq ) 
are equal and cancel each other exactly. As a result the whole functional Su does not give any contribution to l/ri^. 

The same analysic can be applied to Si. We find 



2 

Si[t,t';r;rt] ^ \ j dti j dt2^I{t 



i-t2,r{h)-r{t2)) 
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-I{ti +t2~-t-t', r{ti) - r{t2)) - I(ti +t2-t-t', r{t + t' - h) ^ r{t + t' ~ ia)) 



(69) 



Averaging over diffusive paths r{t) and talcing the limit t — t' —> +00, we observe that is determined by the first 
term under the integral, the other terms give only the time independent (and therefore irrelevant) constant. Thus we 
get 



1 



+ 00 



dt{I{t,r{t)-r{mr 



(70) 



To find the average over the diffusive paths, we introduce the Fourier transform of the function I{t, r) and replace 
(e-^Hr{t)-r(t'))^^ by 6--^"=" I Then we obtain 



+00 



,3-d 



dt 



diod'^k 



-Im 



-An 



cothf— le-*-*-^^'!*! 
\2TJ 



(71) 



As it was already discussed, eq. (^ includes the effect of both electron-electron and electron-phonon interactions. 

To evaluate we use the expression (|^ for the dielectric susceptibility. For typical metallic systems one can 
usually neglect the first term in the expression for e{uj, k) (p^. Then we find 



Im 



-1 



e(w, k) 



Ana 



(tj2 _c2A:2)2 +c4tj2/i:)2 



(72) 



Here c = Wpi(Z?/47r (7)^/2 jg speed of sound in our system. Possible corrections to ( [72|) may become important only 
in special cases of Id and 2d semiconductor systems, where the capacitance may become important. We will come 
back to this point later. 

The expression ( |7^ ) can be significantly simplified for c/D ^k ^l/l. In this limit we obtain 



Im 



i{u:,k) 



Im 



e(cj,fc) 



where 



and 



Im 



Im 



-1 



e(w,fc) 



eph 



e(cj,fc) 



Dck^ 
8a 



Im 



UJ 

Ana ' 



e(cj,fc) 



eph 



{5[lo ~ cfc) - 5{uj + cfc)). 



(73) 



(74) 



(75) 



Phonons with small wave vectors k ^ cjD are strongly damped. For such k we may put Im i^^jjj^^ — 



C. Results 



With the aid of the above results we can now calculate the decoherence time r^. Let us first take into account only 
electron-electron contribution to e and obtain the result for a quasi-one-dimensional system with a ^l. Substituting 
( [7^ ) into (|7^) and integrating over time and the wave vector we arrive at the integral OMr to which diverges at both 
low and high frequencies. The low frequency divergence is cured in a standard mannertj by neglecting the effect of 
environmental fluctuations with frequencies below l/r^. At high frequencies the integral should be cut at the scale of 
the order of the inverse transport time because at higher uj the approximation of electron diffusion becomes incorrect. 
Then we obtain 

, l/re 

1 _e^V2D r duJcoth{Lu/2T) 

T^m ~ J 2^ • ^ ' 
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Eq. (|7|) yields 



, . . , . (77) 

At sufficiently high temperature the first term dominates and the standard result0 ~ (cri/e^D^/^T)^/'^ is recovered. 
This is a classical contribution to r^. As T is lowered the number of classical (low frequency) modes decreases and 

eventually vanishes in the limit T ^ 0. At T < Tq^^ ~ l/^r^Tg the expression ( [tT] ) is dominated by the second term 
and T,p saturates at the value 

w nai/e'^VF (78) 
The estimate for the crossover temperature Tg^^ is obvious from 

-(1 



T^^> ~ evF/2y/ail. (79) 

Making use of eqs. ( |65| , [78| ) it is also easy to find the weak localization correction Sai to the Drude conductance in 
the hmit T = 0. For T < T,^^^ we obtain 

— = - — (80) 

i.e. Sai ~ —cTi/y^Nch, where N^h ^ PpS is the effective number of conducting channels in a Id mesoscopic system. 
For 2d and 3d systems the same analysis yields 



2 

[i + 2rTein(rT^)], 2d, 



1 

— = + Te ln('J „ 



1 



„2 

N3/2 



l + 6iTT,r^% 3d, (81) 



According to (|81|) in 2d and 3d systems the decoherence time becomes independent of T already at relatively high 
temperatures. In the 3d case such temperatures can be of the order of the inverse transport time, while for 2d systems 
we have ^ - (2re ln(p|,aZ))-i. 

Note that in the 2d case we again assumed a < /. Only provided this condition is satisfied the above results for Id 
and 2d systems ave valid for the whole temperature range. At sufficiently high temperatures this condition can be 
softened because in this case r^p is determined by the low frequency fluctuations of the environment. Then the system 
can be considered as a quasi-ld (quasi-2d) one if its transversal dimension is smaller that the corresponding phase 
breaking length a ^ ~ ^JDt^. However at low temperatures high frequency modes become important and the 
situation changes. E.g. in the limit a ^ I the diffusion process has two (and sometimes even three) stages: at short 
times (i.e. at frequencies higher than D/a?) diffusion is obviously 3d, whereas for longer times it can be Id or 2d. 

For a ^ the system is obviously 3d at all T. In the intermediate case I a ^ L^p one should use the 



corresponding Id or 2d formulas for r^p at high temperatures T 3> Tqf^ and the 3d result (81) in the low temperature 
limit T Tqf^ . The crossover temperature Tqf^ can be determined either directly from the integral (^l]) or just by 
comparison of the corresponding expressions for r^. We get Tqf' « Tq'^\a/iTl)'^^'^ . Thus in the limit a 3> I the low 
temperature saturation of r^p takes place already at somewhat higher temperatures than in the case of small a. The 
saturation value of becomes smaller for large a since it is defined by the 3d result (^Tj) . Such details may be quite 
important for a quantitative comparison with experimental results (see below) . 

The above discussion applies to metallic systems. In the case of two-dimensional electron gas in semiconductors 
with I a Lip the same arguments hold, but the saturation value at low T is given by 2d (not 3d) result. The 
classical-to-quantum crossover takes place at T ~ Tq^\a / AttI) . 

Now let us analyze the effect of the electron-phonon interaction on the decoherence time r^p . Substituting ( |75| ) into 
ri|) and applying the same cutoff procedure after a simple integration one obtains 

InU^ +1 T, Id, 



^eph 4t^2^^ V3c/ inaic 
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TT— + 7r5^7^lMmin , 2d, (82) 



eph STrVaWFTe 4712(72 V [Sc' j 



, Tl '. 
T min < 1, — > , 3a. 



Comparing these expressions with the above results for r^p we find /t^^ ~ c/vp ^ 1, i.e. at low T the decoherence 
effect due to zero-point fluctuations of lattice ions is always much smaller than that of quantum fluctuations in the 
electron bath. The temperature dependent term in the expression for the inverse dephasing time due to phonons is 
of the same order as the corresponding electro-electron term in Id and 2d, and can be even bigger in 3d at relatively 
high temperatures T > c/l. Thus we can conclude that in Id and 2d systems the decoherence effect due to phonons is 
never important, while in 3d systems it can modify the temperature dependence of the decoherence time at sufficiently 
high T. 

The above results in Id and 2d are valid if the number of conducting channels Nch in the system is sufficiently big. 
In typical metallic systems this condition is usually well satisfied. However, in semiconductors one can in principle 
achieve the situation with Nch ^ 10. In this case the first term in the expression for e ( ^6|) cannot be neglected in 
general. Moreover for Id and 2d samples with small Nch the energy of the electromagnetic field outside the sample 
may also give a substantial contribution. In order to account for thhis effect we inti-oduce the effective capacitance of 
the system C. Then the influence functional for the field V has the form (see alsdj) 



iSc = ■ 



C(cj,fc) fc2(e(w,fc)-l) 



2 47r 

y 72^^'(-"'-fe) 4. ^ ^"th(-)y-(.,/c); (83) 



y+(w,fc) 



2 



where C{oj,k) ~ (1 -I- es(a;))/41n(l/fca) for a Id wire and C{u!,k) = (1 -I- es{uj))k/8TT for a 2d film. Here es{uj) is the 
dielectric susceptibility of the substrate. The Fourier transform of the function I{t, r) takes the form 



_ /: + us _ WCOth(^) 



7k^ + '^dQ 

Substituting this expression into (fol) we get for Id wire 



3^V2D f doj coth(w/2r) / / /(w) 



,(T) ai J 2^ V^(l + 2/(c^)) V \l + f{to) 



(85) 



where /(w) — 4o^^"^in°(i/fco^) ^^"^ the value koa is roughly of order one. We estimate 

/(.) . i±M^ . {l + e.i.))ppr, ^ 
e^iVos Nch 

where tb = l/me^ ~ 0.5A is the Bohr radius. For metallic wires ppfB ~ 1, Nch S> 1 and the function /(lu) is usually 
small unless £^(0;) ^ 1 at frequencies of the order of I/tc- However for semiconductors f{uj) may be large and r^p 
may become significantly longer than one could expect from eq. ([77|). 

The same analysis can be carried out for 2d films. In this case the effect of capacitance is described by the function 

{1 + es{u}))^/Duj {1 + es{uj))y/uJT^rB 

h[^) = 5 ^ ; ■ 

87rCT2 pfI'O, 

Again one can conclude that this effect is typically negligible for metallic films. For semiconductors with small Nch 
the above effect might cause an increase of r^. 
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V. QUANTUM KINETIC APPROACH AND LANGEVIN EQUATION 



Let us now demonstrate how the usual quantum kinetic description can be derived from our analysis. We start 
from the equation for the density matrix pv (p5|). Rewriting this equation in the "interaction representation", i.e. 
substituting pv e^^^°^ pve^^°* we find 

= -eV^it)pv + PveV+lt) - |((1 - Pv)V-{t)pv + PvV-{m - Pv)), (87) 

where V^{t) — e'^^"*V^ {t)e^^^°* . Let us integrate this equation over time, then substitute the resulting expression 
for Pv into the right hand side of eq. ( p7| ) and average over V^. If the Coulomb interaction is sufficiently weak one 
can proceed perturbatively in V and neglect the dependence of the density matrix pv on this field in the right hand 
side of the resulting equation. Then the result of averaging can be expressed in terms of the correlation functions 
(VV). More precisely, two such functions turn out to be important: 

{V+{h,ri)V+{t2,r2)) - /(ti - h,ri - rs), 

(F+(ii,ri)F-(t2,r2)> =zi?(ti-t2,ri-r2), (88) 

The function {V^{ti,ri)V^ {t2,r2)) differs from zero only for ti > t2- The correlation function (V^V^) is zero for 
all times. Taking this into account we obtain 



dp 

dt 



t 

J dt'(^ V+(t)V+(t')p{t') + V+{t)p{t')V+(t') + V+{t')p{t')V+{t) - pit')V+{t')V+{t) 



-oo 



■ix>+(t)(l - p{t'))V-{t')p{t') - ^V+{t)p{t')V-{t'){l - Pit')) + 1(1 - p{t'))V-{t')p{t')V+{t) 



+lp{nv-{t'){l-p{t'))V+{t)) , (89) 
^ ' v+,v- 



where p{t) = pv=o{t)- 

For simplicity let us consider a clean metal. Making use of the momentum conservation one can significanly simplify 
the equation (|89|). In this case the density matrix is given by p{ri — r2) = / ^J!^3, npe^'P^'^^^'^^\ The operator e"'^"* 
reduces to e^*^''*. Performing the averaging with the aid of eqs. (^,^6|j5^) we find 



dup /■ , _ / — 1 



dt 



dujd klm 



( fc2e(a! fc) ) ^ ^^"'^ ^ (^) ''^^^'^ ^ ~ ^^'^^ " ^^^''^ ^ np^k{l - rip) 



(90) 



The right hand side of this equation represents the standard collision integral which vanishes in equilibrium, i.e. for 

np = l/(exp(ep/r) + 1). 

The equation (^0|) can be also rewritten in the fohowing form: 



^ = ^ / du^Jd^klm (^2^) [si- + ^P-k - ep)iV.np-.(l - Up) - 





~S{uj + Cp-k - Cp)(l + ^i^)"p(l - "-p-fc) + ^{^ + Cp - Cp-fc)(l + N^)np^k{\ - Up) - 

S{uj + £_p - £_p^k)N^np{l - np_k)>, (91) 



where N^^ = 1/ (exp(ti'/T) — 1) is Bose function. This equation describes the standard photon absorbtion and emission 
processes and thus establishes a transparent relation between our approach and one describing the kinetics of an 
electron interacting with the quantized electromagnetic field. In our case the field V is due to fluctuations of conducting 
electrons (or lattice ions - see below). It is quite clear however that the physical nature of this field is not important 
for the electron dynamics, at least as long as this Bose field remains in equilibrium. 

It is important to emphasize that the effect of electron-phonon interaction is also taken into account in eq. (|9l]). 
The phonon spectrum is determined by the equation e(aj, fc) = 0, i.e. the function ^(^\-^ has a pole at a; = ujph{k) ~ iO. 
Therefore one can write 
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Im 



-1 

-1 



A{k) ^ 
= A{k)6{L0 ~ ujph{k)) + 



(92) 



where other contributions to e^^ are denoted by dots. The value A{k) determines the strength of electron-phonon 
interaction. Within the simple model ( [Tq) one has A{k) = Dck^ j'&a — ck^ jX^e^N^. 

Substituting the expression ( |9^ ) into ^^) and integrating over lo we reproduce the standard electron-phonon collision 
integral: 



'eph 



2e2 

r2 



S{ujph{k) + ^p_fe - ^p)iV^^,^(fe)rip_fc(l - rip) 



(5(wp,,(fc) + ^p_fc - $p)(l + 7V„p^(i.))np(l - np_fc) + 6{ujph{k) + £,p - Cp-fe)(l + A''a;p;,(fc))"p-fe(l - "p) 



-(5(wpft(fc) + - Cp-fc)^c^p;.(fe)"-p(l - '^P-fc) ^ 



(93) 



This result demostrates that the function e{uj, k) correctly describes both electron-phonon and electron-electron inter- 
actions. It is not surprizing, because this function just accounts for the collective effect of the environment. Electrons 
propagating in a metal "feel" only the fluctuating electric field produced by the environment, both by electrons and 
lattice ions. Therefore it is quite natural that both contributions can be successfully treated within the same approach. 

The equations (|90|-p3|) are applicable if the distribution functions Hp and Ni^ are close to the equilibrium Fermi and 
Bose functions. It is not difficult to generalize this approach for stronger deviations from equilibrium. Actually the 
electron-phonon collision integral ( p3| ) remains the same in this case, only the distribution function can deviate 
far from the Bose function. In order to generalize the electron-electron collision integral we make use of the following 
nonequilibrium formulas: 



Im 



-1 



e{uj,k)J 7rfc2|e(tj,fc)|2 



fc4|e(a;,fc)|2 



(fp6(uj - £,p+k + ^p){np+k - np), 



d^p[np+k{^ - Up) + np(l - np+k)]5iuj - ^p+fe -|- ^p). 



(94) 



which can be easily derived from (^5|) and ( B6 ) respectively. Substituting these expressions into the equation ( pO[ ) we 
arrive at the electron-electron collision integral for the degenerate plasma: 



f d^k 


d^p' 




1 (2^)3 


(27r)3 


V J 



7 r^TH^i^p'+k + Cp-k - 4p' - Cp) X 

t,p+k — ^p, K)\ 

X [np/+fcnp_fc(l - ?T-p')(l - Hp) - npiUpil - np'+fc)(l - Up^k)] 
Thus the kinetic equation can be written in a standard form 

dur, 



dt 



-^eph ~^ lee ■> 



(95) 



(96) 



where the collision integrals I^ph and I^e are defined respectively by eqs. ( |93| ) and (|95|). 

In order to estimate the characteristic electron scattering time we have to substitute the function Up + Sup instead of 
Hp in the collision integral (|9^). The inverse inelastic scattering time I/ti is then defined as a coefficient in front of the 
term Sup describing deviations from equilibrium. Making use of an obvious identity 1 — 2np_fc = tanh((^p — lu)/2T) 
we get 



Ti [p) TT^ / '^^'^ ( k'^e{u), k) 



5{uj + ip-k - Cp) ( coth — -I- tanh 



2T 



(97) 



It is clear from this equation that the time Ti becomes infinite at zero temperature and at the Fermi energy due to the 
Pauli principle. The same is true for the inelastic scattering time due to electron-phonon interaction. Thus within the 
same general formalism we have demonstrated a clear distinction between the decoherence time and the inelastic 
scattering time r^. The latter becomes infinite at T = because the excitation due to noise (coth(a;/2T)) and the 
energy losses due to radiation (tanh((^ — uj)/2T), see the discussion below) exactly compensate each other. This is in 
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the agreement with the PauH principle: electrons cannot change their energy because at T = all states below the 
Fermi level are occupied. At the same time the value r^p remains finite even at T = because this time is sensitive 
to random noise of the environment only (be it classical or quantum). The effect of such noise leads to decoherence. 
We will return to this discussion below. 

The above kinetic equations were derived for a simple case of a clean system and do not account for the effect of 
elastic scattering. In the case of a disordered metal the electron mcwientum is not conserved and the whole derivation 
becomes much more complicated. One can demonstrate (see e.g.t3) that in the diffusive limit the result is roughly 
equivalent to a substitution 

(5(w + Cp_fe -Cp) -> Re 



+ Dk'^ 

in the expression (p7[). An extended analysis of the inelastic scattering time in various limits is given in Rcf. 

One can also formulate an alternative approach and derive the quasiclassical Langevin equations describing electron 
dynamics in a disordered metal. In doing so, we follow the same procedure as one described in Refs. |5|,p9|. 

Consider only close electron paths for which the values r — ri — r2 and p" — pi — P2 are small. Then we can 
expand the effective action in the exponent of eq. (^3|) in powers of r~ and keeping only the quadratic terms. 
The action becomes Gaussian in terms of these variables and the integral ( p3| ) is dominated by the saddle point 
trajectories: SS/Sp" = and SS/6r^ = 0. The first equation coincides with one without dissipation: r — p/m. With 
the aid of this equation the momentum can be easily excluded and we get: 

t 

mr + \/U{r) + e^[l~2n{r,mr)] J dt'\/rR{t ~ t' ,r{t) - r{t' j) ^ ~eE{t,r). (98) 

— OO 

Here E{t, r) is a fluctuating electric field. Equilibrium fluctuations of this field are described by the correlator 

{m„r,)E,it,,r2))=A.S,,J ^|gp:inr(^^^^)e--(*-*^)+»'^(-~-) (99) 



If needed, the generalization of (98|99|) to a strongly nonequilibrium situation can be also provided. Also more general 



expressions for the kernel i?(i, ri,r2) and for the correlation function (M) for the case |ri — < I can be easily 



derived. Combining these expressions with (98) one can obtain a quasiclassical description of electron dynamics also 
at scales ^ /. However such details are not important for us here, eq. ( p8| ) is presented merely to illustrate important 
physical effects. 

The equation (^ obviously satisfies the requirement of causality and captures all the essential features of electron 
dynamics in a metal. E.g. it demonstrates that electrons in a metal cannot infinitely decrease their e nerg y: effective 
damping due to the presence of the environment (described by the last term in the left hand side of {p8\)) is zero at 
the Fermi energy (n — 1/2) and becomes negative below this energy. Thus electrons with the initial energy above fi 
will loose it before they reach the Fermi level. On the contrary, holes with the initial energy below fi will be pushed 
up to the Fermi surface. This simple example demonstrates again that our analysis accounts explicitely for the Pauli 
principle. The corresponding information is contained in the influence functional which depends on the occupation 
numbers. 

The damping term in eq. ( |9^ ) depends on the function R{t, r) which is determined by the correlation function 
{V^V~) (see (p^). The physical origin of this damping term is quite transparent: the electron (or the hole) prop- 
agating in a metal produces the screened electric potential due to the presence of other electrons and "feels" this 
potential itself. In this sence eq. ( |98| ) is similar to the equation of motion of a high energy particle (e.g. muon) in a 
metal. The important difference between these two cases, however, lies in the factor 1 — 2n which is present in our 
case due to the Pauli principle. Formally this factor enters due to fluctuations of the field V~ which is "sensitive" to 
the Pauli principle. The fluctuating electric flcld E in the right hand side of eq. (Bq) is, on the contrary, not affected 
by the Pauli principle because its correlation function depends only on the field V^. 



With the aid of eq. (M) it is also easy to understand why the real part of the influence functional Sa (p4) does 



not contribute to the decoherence time. According to (98) the phase difference acquired by the electron propagating 



along some classical path can be split into two parts: the regular contribution due to damping {Sr) which depends 
only on the electron trajectory, and irregular part due to noise (Si). Considering now the contribution from a pair of 
time reversed paths, we observe that the regular contributions are the same and cancel each other because they enter 
with a different sign. Only irregular contributions due to noise survive and determine Tip. For each path the regular 
contribution may have a different value depending on the path and energy and even vanish (for energies at the Fermi 
level). However by no means this affects the noise terms and thus which always remains finite. 
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VI. DISCUSSION 



With the aid of a general formahsm of Green-Keldysh functions we have analyzed a fundamental effect of quantum 
decoherence of the electron wave function in a disordered metal due to interaction. Our treatment was carried out with 
no more assumptions than the usual ones in the weak localization theory: the elastic mean free path was considered 
large as compared to the Fermi wavelength ppl ^ 1 and interactions were assumed to be sufficiently weak. Actually 
the main effect does not really depend even on these general assumptions, it can be observed already from a formally 
exact equation for the density matrix ( ]2^ , p5| ). 

We have demonstrated that the effect of interaction of the electron with other electrons and lattice ions in a 
disordered metal is equivalent to the effect of a dissipative environment. This effect is similar (although not exactly 
equivalent) to that of a dissipative Caldepa-Leggett bath. Fluctuations in the dissipative environment play the key 
role in the effect of quantum decoherenceES. 

Although the environment is dissipative it is not inelastic scattering processes with energy transfer between real 
quantum states that cause quantum decoherence at low T. Rather this is the effect of quantum noise in a dissipative 
environment. Due to interaction with quantum fluctuations of this environment the electron not only "goes" through 
virtual states with different energies but after a finite time also looses information about its initial phase. At the 
same time no real processes with energy transfer exist at T = 0. ft is also interesting to note that - as can be seen 
e.g. from eq. ( |8^ ) - voltage fluctuations and therefore the decoherence rate l/r^ vanish in both limits a ^ and 
cr — *■ OO. 

The presence of Feynman paths which strongly deviate from the classical paths depicted in Fig. la is also crucially 
important for the whole effect. In the case of disordered metsLls-Lhese are time reversed paths (Fig. lb). In the case 
of Caldeira-Leggett models other paths (e.g. instantons, seeE30'El) play a similar role. If the important paths go 
sufficiently far from each other, electrons gain an additional phase due to the environmental noise. This phase always 
vanishes for the classical paths of Fig. la. 

We clarified the role of the Pauli principle in the effect of quantum decoherence. Although the Pauli principle plays 
a crucial role in the inelastic electron-electron collisions it does not affect quantum decoherence. The physical reason 
for that is transparent: fluctuations of the collective electric field in a metal are equivalent to local fluctuations of 
the Fermi energy or the electron density (see Section 2). Such fluctuations are essentially the same for many electron 
states, and the Pauli principle does not play any role. 

On a very simple level the decoherence effect at low T can be understood as follows. During interaction with other 
electrons the electron phase changes. These "other" electrons (or some of them) "go out of the game" (the system is 
open!) and after some time the information about the initial electron phase is lost. The same mechanism works for a 
particle interacting with the Caldeira-Leggett bath. 

We would like to emphasize that no special assumptions about the "open system" were made in our analysis. We 
just used the standard formalism of Green functions applicable to any grand canonical ensemble in condensed matter 
physics. Already by fixing the chemical potential and allowing the total number of particles to vary one defines 
the system as open. Thus our results are general and in principle can be applied to any disordered system perhaps 
except for ultrasmall objects with fixed numbers of particles. The latter appear to be irrelevant for weak localization 
measurements. 

Our analysis provides no support for the point of view according to which the effect of high frequency (quantum) 
fluctuations of the electronic environment on should be detemined by inelastic processes with high energy transfers 
and therefore should be obtained from the kinetic analysis. The kinetic treatment deals with paths of Fig. la and 
allows to calculate only but not r^. In contrast, the latter is determined by time reversed paths of Fig. Ic for 
all relevant frequencies. In this sence there is no qualitative difference between the effect of classical and quantum 
fluctuations of the environment on r^. Both should be treated on equal footing. 

There is also no evidence that the Pauli principle "forbids" decoherence at low T. Therefore the suggestion to 
account for the Pauli principle by introducing an additional term proportional to tanh(ct;/2r) into the expression for 
Tip does not appear to be justified. We have demonstrated that the combination "coth — tanh" enters the expression 
for 1/Ti (97) but not for the inverse decoherence time 1/r^. The latter contains only "coth" and does not vanish at 

We believe that the same should hold for a sometimes conjectured cancellation of Keldysh diagrams at T = 0. 
To the best of our knowledge a rigorous formulation of the diagramatic representation ictr is still lacking in the 
literature. An attempt to provide this formulation was made by Fukuyama and AbrahamsEj. But it ia-well understood 
by nowEa that only diagrams relevant for the inverse inelastic time 1 /tj were taken into account inEj. Cancellation 
of these diagrams at T = is quite natural but it does not yet tell anything about . Since our analysis should in 
principle include all diagrams we see no reason to expect any cancellation of diagrams for even at T = 0. 

The low temperature saturation of can cause dramatic consequences for the existing picture of strong 
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localization0ii in lower dimensions. There seems to be no reason to challenge this picture as long as one neglects 
Coulomb interaction between electrons. But - as was demonstrated above - it is exactly this interaction that leads to 
quantum decoherence at T = 0. If one compares the effective decoherence length L^p with the localization length Lioc 
one immediately sees that for the parameters of a typical metal the former is always smaller than the latter. E.g. in 
Id our analysis yields L^p ^ yj Dr^p ^ l\/ Nch- This expression is parametrically smaller than the localization length 
Lioc ~ INch for large number of conducting channels. Analogously in 2d at low T the weak localization correction to 
conductivity saturates at the level 5a2 ~ — (e^/27r^) ln(p|,aZ) and the crossover from weak to strong localization never 
occurs. In other words, all that implies that (at least for pp^ ^ 1) due to electron-electron interaction the electron 
will loose coherence already before it can get localized, i.e. strong localization does not take place at all and the Id 
and 2d metals do not become insulators even at T = 0. 

We can also add that it would be interesting to investigate the effect of electron-electron-jrteraction on more 
sophisticated weak localization corrections, like ones due to quantum interference of diffuson£3. Here the point is 
that (unlike Cooperons) diffusons do not decay at T = even in the presence of interaction. Nevertheless, we cannot 
exclude that even at T = quantum interference of diffusons can be also destroyed due to the same mechanism of 
quantum noise in a dissipative electronic environment as it was discussed above for Cooperons. This problem will be 
studied elsewhere. 

Irrespectively to such details it is quite clear that an adequate description of electron-electron interaction can 
seriously change the existing picture of strong localization in low dimensional disordered metals. This perspective 
is very interesting but not very suijprizing. In fact, the important role of electron-electron interaction was already 
pointed out long ago by FinkelsteinEa. 

One more comment concerns the relation between the low temperature saturation of the decoherence time r^p and 
the applicability of the Fermi liquid theory. It is sometimes believed that the applicabitity condition for the latter is 
Ttcp ^i-il- From a practical point of view this condition is not very restrictive: e.g. for typical experimental values 
of Ref.El the value was found to be in the range 10~^ 10 ns and the above condition can be violated only at 
temperatures below 1 100 mK. However theoretically it is very important to understand if a finite decoherence time 
at T — > would imply the breakdown of the Fermi liquid theory at low T. We believe that there is no direct relation 
between Tip and the lifetime of quasiparticles. The latter should be rather defined by some inelastic time (e.g. by r^) 
which is always much longer than at low T. Thus it appears that even for Tt^ < 1 there is still no evidence for 
the breakdown of the Fermi liquid theory in metals. Previously an analogous conclusion was reached in Ref. 

Finally, we briefly discuss the agreement between our results and the available experimental data. The comparison 
between theoretical and experimental values for the decoherence length L^p = ^/Dt^ at zero temperature is given in 
the table |. To calculate L^p we first estimate the decorence time with the aid of ([77|). At T = the time r^p can be 
conveniently expressed in terms of measurable quantities: 

r^p = V6|^, (100) 

where Rq = 7r/2e^ = 6453 is a quantum resistance, L is the wire length and R is the total resistance. The Fermi 
velocity for gold wires was taken to becil 1.39 x 10^ m/s. The diffusion coefficient D was estimated with the aid of 
the Drude formula a = 2e^A^o-D. The density of states fo*; gold is chosen to becJ A'^o ~ 6 x 10^^ s/m^. Note that the 
numerical values for D are not identical to those given iru. 

M 



TABLE I. Theoretical values of the coherence length = -y/ Dt^ in comparison with the experimental results 



Sample 


w, ran 


t, ran 


L, lira 


R/L, n/fj.m 


D, 10"^ mVs 




rthoor 


Au-1 


60 


25 


57.9 


29.14 


7.8 


5.54 


1.8 


Au-2 


110 


60 


207 


1.46 


35.5 


16 


16.5 


Au-3 


100 


35 


155 


9.31 


10.5 


5.2 


3.6 


Au-4 


60 


25 


57.9 


31.29 


7.3 


3.6 


1.6 


Au-5 


190 


40 


18.9 


191.7 


0.24 


0.35 


0.12 


Au-6 


180 


40 


155 


2.91 


16.3 


8 


8.1 



The width, the thickness and the length of the wire are denoted respectively by w, t and L; R/L is the resistance 
per unit length and L^^, L^^^°'^ are experimental and theoretical values of the decoherence length. The agreement 
between both looks very reasonable for all samples, especially if one takes into account an uncertainty in a numerical 
prefactor in our formulas due to the cutoff procedure and possible effects of the sample geometry. It is also important 
to emphasize that our comparison involves no fitting parameters. 

The temperature dependence of the inverse decoherence time l/r^p is plotted in Fig. ^ The agreement between 
theoretical and experimental results is again reasonable. It appears, that within the temperature interval of Fig. 
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temperature dependence of the experimental data is not far from a linear one. Note, that the temperature dependent 
part of eq. is sensitive to details of the low frequency cutoff. The plotted theoretical dependence corresponds to 
the lower cutoff frequency equal to \|T^p at T = 0. With a slight adjustment of this cutoff one can in principle reach 
a perfect agreement with the experimental data. However, we believe that this adjustment is not needed: within the 
accuracy of our calculation the agreement is already very good. 




500 



2500 



T, mK 

FIG. 2. The dependence of the decoherence rate 1/t<^ on temperature for one of the thin gold wiresEJ. The wire has the 
following parameters: w — ?>b nm, t = 20 nm, L = 231 /im, R = 4411 f2, D = 2.56 x 10~^ m^/s. Experimental points are 
shown by squares, the solid curve represents the teoretical dependence ([771). The experimental data were kindly provided by 
P. Mohanty. 



Although we do not present the comparison of our theoretical results with the experimental data from Refs. 21-23 we 
note that all these dataj-are sufficiently well described by our formulas ( [77|j8l| ). The agreement with the experimental 
results by Mueller et a/.EJ for 2d gold films turns out to be worse: the value estimated from eq. ( pT[ ) is considerably 
shorter than the experimentally measured value. The reason for this discrepancy is not clear at the moment. We 
can speculate that it might be due to the high purity j«f the filmsEs which can make the diffusive approximation 
insufficient. If this is the case the resistance of the filmsEil is rather due to electron scattering at the boundaries and 
our theoretical analysis should be substantially modified. This problem will not be discussed here. 



pm weak to strong localization was not detected even 
I a sharp increase of the resistance by several orders 



In all the above (and many other) experiments no crossover f 
at the lowest temperatures 20 ~ 40 mK. Very recently in Ref.t 
of magnitude was observed in quasi-ld GaAs structures with small number of conducting channels Nch ^ 10. This 
increase of the resistance was interpreted as a crossover from weak to strong localization. Although this could be one 
possible scenario we believe the interpretation of the above effect is still an open problem. Indeed, if one estimates 
the crossover temperature from the standard criterion L^p Lioc and makes use of (|l|) for the parametergla one 
arrives at the crossover temperature T ~ 1 10 mK. However the increase of the resistance was observed at a much 
higher temperature T > 1 K. This discrepancy might indicate the important role of Coulomb effects in the system. 
Analogously Coulomb effects were found to play the key role in many experimens with disordered 2d systems, see 
e.g.E3 for the corresponding discussion. Further experiments are definitely needed for better understanding of the 
electron-electron interaction effects in disordered low dimensional systems. 

We would like to thank Ya. Blanter, G. Blatter, C. Bruder, D. Geshkenbein, V. Kravtsov, A. Mirlin, A. van Otterlo, 
M. Paalanen, A. Schmid, CpSchon and P. Wolfle for valuable discussions, comments and/or encouragement. We are 
also indebted to the authorso for communication of their new results prior to publication. This work was supported 
by the Deutsche Forschungsgemeinschaft within SFB 195 and by the INTAS-RFBR Grant No. 95-1305. 
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APPENDIX A: DENSITY MATRIX 



With the aid of eq. (|9|) the time integral in the last term of the Dyson equation ( |l8| ) can be transformed as follows: 
t t 
dt'Go{ti,t')eV{t')G{t', t2) = G^ihM) - J dt'G„{h,t') (^t-^ - iJo(i')) ia) 



= Goih,t2) - tGo{h,t)G{t, h) + 2Go(ii, 0)G(0, ta) + 
dt' O-^+Hoit')) Go{tut')G{t',h) 



= GoihM) - G(hM) ' tGo{ti,t)G{t, h) + iGo(ti, 0)G(0, ta) 
Here we performed the integration by parts over the time t' and made use of the equation 

d 



(Al) 



which defines the field- free Green-Keldysh function Go. Substituting the result ( |Al[ ) into eq. ( p^ ) we arrive at (|2 
Let us substitute the representation into (^2|). Then we find 



Mh,t)i-b + h{t))Uv{t, t2){a + fv{t2)) - MhM^ + /o(0))C/y (0, h){~b + fv(t2)) = 
It is now convenient to represent the operator fv in the form 

fv{t2) = Uv{t2,0)git2)UviO,t2). 

Then for g we get 

[l + {b- fom{Uo{0,t)Uvit,0) - imt2) = /o(0) - (S- /o(0))({7o(0,<)C/y(<,0) - 



(A2) 
(A3) 
(A4) 



Let us put t2 = t and introduce the scattering matrix S — ?7o(0, t)Uv{t, 0) = sia + S2b. This matrix is diagonal because 
both C/q a-iid Uy are the diagonal matrices. The matrix elements si and S2 are defined as si^2 — uo(Oi ^)wi,2(i, 0). 
Making use of the above notations and rewriting the equation (A4) in components we find 



1 
1 



-l + Po 1-po 



si-l 

S2 - 1 



1 - Po Po 



Po -Po 
-1 + po 1-po 



Si - 1 
S2 - 1 



911 912 

921 922 

I 







(A5) 



Multiplying matrices and keeping only the part of the resulting matrix equation which depends on 1712 and 322 we 
obtain 



l+Po(si-l) -P0(S2 - 1) 

-Si + I + po(si - 1) S2 - Po(s2 - 1) 



512 
.922 



(A6) 



Substracting the second equation from the first one we get Sigi2 — 82922 — 0, or 522 — ^2 ^^1(712. Subsituting this 
result into the first equation we find 



[1 +P0(S2 ^Si - l)]5i2 = 



Pa- 



(A7) 



Note that the S-matrices enter the above equation only in the combination S2 "'^si which does not contain the evolution 
operator uq, i.e. we have s^^si = U2{0,t)ui{t,0) and 



[1 - po(0)(l - U2{0, t)uiit, 0))].gi2(t) = po(0). 



(A8) 



Rewriting (A3) for the 1,2-component of the matrix fv{t) 

fl2{t)=Ui{t,0)gi2{t)u2{0,t), 

and making use of the identity py(i) — fi2{t) we arrive at the result 
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APPENDIX B: EFFECTIVE ACTION AND FDT 



For Vi_2 the electron Green functions G12 and G21 can be expressed in the form 

G12 = iuo(ti,0)po"o(0,<2), G21 = -mo(ii,0)(l - /9o)wo(0, <2)- 



(Bl) 



In thermodynamic equilibrium we have uo(i, 0)^0^^0(01 1) — Po for any time t. 

Let us introduce the basis of the eigenfunctions for the single electron Hamiltonian, Hoijjk — ^kipk- Without loss of 
generality we can choose these eigenfunctions to be real. The initial density matrix po is assumed to be diagonal in 
the basis ipk, namely po = X^fc "-fclV'fc) (V'fcl- This assumption is justified only for weakly interacting particles. Then 
the functions G12, G21 (Bl), G^ and G^ (29) can be written in the form 



Gi2{ti,t2,ri,r2) = i^e-*«'=(*i-*^)nfcVfc(r-i)V'fe(r2), 
fc 

G2i{ti,t2, ri, ra) = ^ e-*«^(ti-*2) „ nk)Mri)Mr2), 

k 

G«(ti, i2, r,,r2) = -^0(^l - ^2) ^ e-^^''^*^"*^ Vfc(r-i)^fc(r2), 

k 

G^{h,t2,ri,r2) - i0{t2 - h) J2 c~'^''^''-''^Mri)Mr2)- 



With the aid of these expressions the kernels (gljJS^) can be represented as follows 

X{t, ri, ra) = 2te^9{t) ^ e-*(«'--«')*(nfc - ng)Mri)Mr2)Mr2)Mriy, 

k,q 

r]{t,ri,r2) = ^e^^e-'^^''-^^'>*[nk{l - Ug) + ng{l - nk)]'^k{ri)iJk{r2)^q{r2)'tpq{ri). 

k,q 



Performing the Fourier transformation in time: 
X(^,n,r,)^-2e^^ 



k,q 



UJ-£,k+£.q + iO 



TT:ipkiri)ipk{r2)i'qir2)tpq{ri), 



(B2) 



(B3) 



(B4) 



(B5) 



r]{uj,ri,r2) = 7re^^(5(w - £_k +C«)K-(1 - ng) + n,(l - nk)]tpk{ri)i^k{r2)-ipq{r2)'ipq{ri), (B6) 

k,q 



and substituting the equilibrium distribution function rik — l/{e^''^'^ + 1) one immediately arrives at (t 
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